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Coupled phase-oscillators are important models related to synchronization. Recently, Ott- 
Antonsen(OA) ansatz is developed and used to get low-dimensional collective behaviors in coupled 
oscillator systems. In this paper, we develop a simple and concise approach based on the equations 
of order parameters, namely, order parameter analysis, with which we point out that the OA ansatz 
is rooted in the dynamical symmetry of the order parameters. With our approach the scope of the 
OA ansatz is identified as two conditions, i.e., infinite size of the system and only three nonzero 
Fourier coefficients of the coupling function. Coinciding with each of the conditions, a distinctive 
system out of the scope is taken into account and discussed with the order parameter analysis. 

Two approximation methods are introduced respectively, namely the ensemble approach and the 
dominating-term assumption. 


I. INTRODUCTION 

Understanding the intrinsic mechanism of collective 
behaviors of coupled units has become a focus for a vari¬ 
ety of fields, such as biological neurons, circadian rhythm, 
chemically reacting cells, and even social systems hhzi- 
Some properties of collective behaviors depend on the 
complexity of the system, while the other properties, such 
as phase transition, may be described by low-dimensional 
dynamics with macroscopic variables mum- Discov¬ 
ering the method to simplify the system is just as impor¬ 
tant and as fascinating as the discovery of the complexity 
of it. 

Like most cases in physics, simplification and low¬ 
dimensional reduction are associated with some symme¬ 
try of the system. As the identity of gas particles is the 
foundation of statistical mechanics and collective vari¬ 
ables as temperature and pressure, the identity of the 
coupled units in a complex system is also related to some 
order parameters. In previous works, in the limit of 
large number of oscillators TV —* oo with special coupling 
function this work has been done by Ott-Antonsen(OA) 
ansatz [9] for oscillators with nonidentical parameters. 
As for the oscillators with identical parameters the same 
result was also got from group theory analysis in m 
called Watanabe-Strogatz’s approach. 

In this paper, we will show a different way in which 
the low-dimensional reduction is a natural consequence 
of the symmetry of the system by taking the order param¬ 
eters as collective variables, namely the order parameter 
analysis. Our approach is simple and concise for oscil¬ 
lators with both identical and nonidentical parameters. 
We can also get the scope of the OA ansatz. Two more 
cases beyond the scope are further discussed with appro¬ 
priate approximations. With our approach we show that 
the OA ansatz can be used beyond its scope with the 
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approximations works. 

The model discussed in this paper is the all-connected 
coupled phase-oscillators. In the first section, the dynam¬ 
ical equations for order parameters are derived. The OA 
ansatz is got naturally from the symmetry of the order 
parameter equations, with its scope as the limit of in¬ 
finitely many oscillators and the condition that only three 
Fourier coefficients of the coupling function are nonzero. 
In the second and third sections, we will consider two 
approximate use of this ansatz beyond its scope, the case 
of a finite-size system and the case of coupled oscillator 
systems with more complicated coupling functions, with 
approximations respectively, i.e., the ensemble approach 
and the dominating-term assumption. 

II. EQUATIONS FOR ORDER PARAMETERS 

The famous Kuramoto model for the process of syn¬ 
chronization attracts much attentions upon it is proposed 
and has been developed for decades. This model consists 
of a population of N coupled phase oscillators {ifij} with 
natural frequencies {uij } , and the dynamics are governed 
by 

N 

ifj(t) =LOj~Y_1 Ajk sin(^ fe -(Pj), j = 1,... ,N. (1) 
k -1 

With the mean-field coupling Ajk = K/N and the defi¬ 
nition of order parameter 

1 N 

“=jvE e, '-> (2) 

3 =1 

the equations Eq. Q can be rewritten as 

=Uj - Yi^ ae ~ 1Vi j = 1,... ,1V, (3) 

where i is the imaginary unit, a is the complex conjugate 
of a. Apart form parameters I\ and ujj, the dynamics of 
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each phase variable pj depends only on itself and a. It 
is the important character of this mean-field model, and 
the order parameter a is always used to describe the state 
of the system, as the system is in synchronous states if 
and only if |a| = 1. 

A more general form of the mean-field model with 
phase oscillators can be written as 

Vj(t) = F (<x,ipj,(3, 7j), j = l,...,N, (4) 


Moreover, for the coupling function F(a,p,/3) is al¬ 
ways 27r-periodic for p, we have the Fourier expansion of 
F(ot, p, f3), 

OO 

F(a,'p,0)= J2 (9) 

j=— oo 

With the the definition Eq. ([8]) , the dynamics of order 
parameters can be got as 


where F(a,pj, (3,~/j) is any smooth, real, 27r-periodic 
function for pj. a = (...«_!, cxq, is the order pa¬ 

rameter with the nth order parameter a n defined as 

1 N 

°“ = jv (5) 

i =i 

(3 = (ySi, fc—fis) is the identical parameter such as the 
coupling strength K which is identical for all the oscilla¬ 
tors. 7 j = is the nonidentical parameter 

such as the natural-frequency {ujj } which is nonidentical 
for the oscillators and usually has a distribution among 
the system. 

Almost all the mean-field models based on the Ku- 
ramoto model belong to this general category. In the fol¬ 
lowing, we will build our approach for this general model 
to explore the conditions in which we can get the low¬ 
dimensional description of the system Eq. ([4]). 


A. Identical oscillators 

To begin with, let us consider the simple case of iden¬ 
tical oscillators, e.g., uij = w* for j = 1, 2,..., N in the 
Kuramoto model, which reads 

= F{a,p j ,0),j = 1,...,N. (6) 

In the limit of infinitely many oscillators N -A oo, let 
pit, p)dp denote the fraction of oscillators that lie be¬ 
tween p and p + dtp. Because each oscillator in Eq.([6]) 
moves with the angular velocity Vj = F{a,pj, /3), the 
single oscillator density p{t, p) obeys the continuity equa¬ 
tion as 

dp d{pv) = 

dt dp ’ (7) 

v = F(a,<p,/3). 

If the phase density p(t, <p) is known, all of macroscopic 
properties of the system can be got through some statis¬ 
tical average, such as order parameters a n as 

p{t,p)e mip dp,n = 1 , 2 ,.... (8) 

If we are only concerned with the collective or macro¬ 
scopic state of the system, as in this paper, the macro¬ 
scopic description of the phase density is equivalent to the 
microscopic description of phases of all the oscillators. 



/*27T 

a n = ine 171 ^ ipchp 

Jo 

p2ir 

= / ine mip F(ct,(p, (3)dip, 

Jo 


( 10 ) 


where n = 1,2,.... Substituting the expansion of 
F(a,ip, (3) into these equations, we have the closed equa¬ 
tions for the dynamics of order parameters as 

OO 

a n = m ^ fj{a,0)aj +n ,n= 1,2,.... (11) 

j=— oo 

On the other hand, from Eq. ([8]) the order parameters 
are exactly Fourier components of the phase density. If 
all the order parameters are known, we have 


P(t, <f) = 7T- 


OO 


1 

27r ^ 

j—— oo 


CKo-e 




( 12 ) 


Hence if we know all the initial values of order param¬ 
eters, with the dynamical equations Eq. (Ill and the 


Fourier transformation Eq. (12), the system is identi¬ 


fied, which is the same as the dynamical equations Eq. 
([6]) for phase variables {pj} and the continuity equation 
Eq. 0 for the phase density p(t,p). Therefore, we can 
choose either the order parameters, the phase variables or 
the phase density to perform the analysis of the collective 
behaviors of the coupled phase oscillators. 

These three descriptions of the system, i.e. phase vari¬ 
ables, density of phase and order parameters, correspond 
to the dynamical, statistical and macroscopic descrip¬ 
tions of the system, respectively, and they are equivalent 
to each other in the limit of infinitely many oscillators 
N —> oo. As a matter of fact, the Ott-Antonsen ansatz 
[9] is based on the representation of the density of phase, 
and the Watanabe-Strogatz’s approach m is based on 
the scenario of the phase variables. In the following, we 
will take our approach on the base of order parameters. 

The complexity of dynamical equations for order pa¬ 
rameters Eq. (Ill depends on the coupling function, 


or explicitly, on the Fourier expansion of F. First, let 
us consider the simplest case, with only the first three 
nonzero terms of the Fourier expansion, 

F{a,p) = fi(a)e 1<p + /_ i{a)e~ llp + / 0 (a). (13) 

Because F is a real function, we have /_i = /i where /i 
means the complex conjugate of /j. With Eq. (13), the 
dynamics of order parameters Eq. 0) becomes 


fl&n—l H - fo&n)i 


(14) 







3 


with n > 0 and a_ n = a n . The recursion form of these 
equations shows that there is some structure for order pa¬ 
rameters with which we can simplify the system and get 
some low-dimensional equations for the high-dimensional 
coupled oscillator dynamics of the system. In fact, by 
choosing an invariant manifold as a„ = a",n > 0, all 
the equations for a n is reduced to 

di = i(/i(ai )oq + /i(ai) + /o(ai)ai) (15) 


for n = 1,2,.... The infinitely many dynamical equa¬ 
tions of order parameters is thus reduced to a single equa¬ 
tion on this manifold and the corresponding phase den¬ 
sity defined by Eq. (12) is the so-called Poisson kernel 
distribution as 


p(t, ip) = 


1 


1 - r 2 


2-7T 1 — 2 r cos (ip — </>)+ r 2 ’ 


(16) 


where op = re. If we choose the initial state on this 
manifold, the state will never evolve out of it, which is 
called the invariant manifold of the dynamical system 
Eq. (14). This is exactly the low-dimensional behavior 


of the system we are looking for. With our approach, 
the manifold is derived naturally and concisely. We will 
call this manifold the Poisson manifold in this paper as 
in m where the dynamics of a\ 1 as Eq. ( |15| , is got from 
the group theory analysis for Joseplison junction arrays. 


B. Nonidentical parameters 

Furthermore, the low-dimensional behavior is not con¬ 
fined to the system of identical oscillators. The order 
parameter analysis we used above can also be used in a 
more general case of nonidentical oscillators with non¬ 
identical parameters. 

Firstly, let us consider a system of nonidentical os¬ 
cillators with a discrete distribution of the parameter 
u> = {a/ 1 ),...a/ p )}. This distribution naturally sep¬ 
arates the oscillators into p groups, and each group has 

the same parameter u)( p \ In this case, the phase variable 

(V) 

for oscillators is denoted by </?)■ which means the jth 
oscillator with the same parameter 

For each group we can define a local order parameter 
as 

«i r) = -E^‘- n=1 ' 2 --- a?) 

n U 

where n r is the number of oscillators in the group with 
oj( r \ The order parameter of the system denoted by a n 
is defined by 


a n='^2a^ ) — ,n = 1,2,.... (18) 

z ' n 

r= 1 

Assume that the coupling function for each oscillator 
depends only on the order parameters whether for 


the whole system or the groups, then equations for phase 
oscillators are 

< pf\t) = F{a,a (r '> 1 (p ( J r) ,uj ir \f3),j = l,...,n r , (19) 

where r = 1,2 ...p. = (...a^, a^, a^...) and 

a = (...at-i, cuo; op...) are the local and global order pa¬ 
rameters, respectively. 

As the coupling function F(ot, /3, a.( r \uj( r \ is 27r- 
periodic for <pj , by performing Fourier expansion of F, 
the equations for order parameters of each group read 

OO 

a ( n r) =m ( 20 ) 

j = ~ oo 


where r = l,2...p, 

E oo 

j — — C 


= 1 , 2 , 




(r) 


and F(ip^) = 


The equations Eq. (20) are closed for local order pa¬ 
rameters, and the local order parameters here can be re¬ 
garded as the coordinates of the system which are equiv¬ 
alent to phase variables. Each group is almost separated 
from others and the only relation for the groups is the 
dependence of the coupling function on order parameters 
of the system. 

In the simplest case when only the first three terms of 
Fourier expansion are nonzero, we have 


a. 


(r) _ 


= + /l04-l + /o«D 


(r) 


vW) 


( 21 ) 


with r = 1,2 ...p and n = 1,2,.... For each group, 
(V) / (r)\ 

an = (oq ) n is obviously a solution for the equations 

(r) 

which reduce all the equations for an J to the same one 
as 


a 


(r) 


= i(/ 1 (aW) 2 + / 1 +/ 0 an, 


M\ 


( 22 ) 


where r = 1,2 .. .p. The Poisson manifold is an invariant 
manifold for each group, and the order parameters of the 
system read 


r—1 


= 5 >i r) ) B ^>»= . 


(23) 


which is determined by a\ r . r = 1,2 ...p. There¬ 
fore, the Fourier coefficients of the coupling function 
depends on a~{\r = 1,2 ...p only. 
The system described by the local order parameters is 
governed by a group of Poisson manifolds with the low- 
dimensional dynamics Eq. (22). 

In the limit with n r —> oo,r = 1,2 ...p and p —> oo, 
the distribution of parameter oj as {77. 1 , n 2 , ...n p } for 
{a/ 1 ), uj( 2 \ ...u/ p )} has the continuous form denoted by 
the function g(ui) and the local order parameters a^\ r = 
1,2.. .p become the function of w, i.e., aq(w). Replac¬ 
ing the summation by the integral over w, we get the 
continuous form of the dynamical equations as 

di(w) =i(/i(w,ai(w),a:*)(ai(w)) 2 + /i(w, Qq(w), a*) 

+ /o(w,«i(w),a*)ai(w)), 


( 24 ) 
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where a:* = aj}, a\...) are functionals of ot\(ui) as 


a * n = J ( a i( w ))" 5 (w)dw, n= 1 , 2 ,.... 


(25) 


Even though for each specific w the approach has al¬ 
ready reduced the dynamics of the oscillators to the dy¬ 
namic of ai(w) as Eq. (24), it is still hard to get any an¬ 
alytical results as Eq. (24) depends on the integral Eq. 


(25) which cannot be expressed by functions of ai(w). 


On the other hand, for some specific choice of gioj), the 
integral Eq. (251 can be obtained analytically, in which 


case we can get the simpler form of Eq. (24). 

What we are looking for is a two-step reduction. The 
first step is finished by introducing the Poisson manifold 
with which we have reduced the dynamic of phase oscil¬ 
lators to the dynamic of a group of Poisson manifolds, 
i.e., Eq. (241. The second step is rooted in the relation 


between the Poisson manifolds in the group which will 
reduce the dynamics of the group further to the dynamic 
of a specific Poisson manifold in this group. 

For instance, in the case that the function aq(u;) is 
analytical which can be extended to the complex plane, 
and the distribution of ui is the Lorentzain distribution 
as 


9{u) = 


7T(( W - W0 ) 2 + A 2)’ 


(26) 


the integral Eq. (|25j) can be obtained by residue theorem 
as a* = (a n (—i)) n . Setting w = —i, Eq. |24j) becomes 


di(-i) =*(/i(—i, «i(—i))(ai(—i )) 2 +/i(-i, aq(-i)) 

+ /o(-i,aq(—i))ai(—i)), 

(27) 

which is closed for order parameter oq(—i) = a\. Then, 
the behavior of a*, together with the collective behaviors 
of the system, is determined by Eq. (27). 


Another solvable example is the case when the distri¬ 
bution of natural frequencies of oscillators is the Dirac’s 
delta distribution function g(tu) = 5(u> — wo). The in¬ 
tegral Eq. (25) can be worked out as a* = (ai(w 0 )) ra . 
Setting to = wo in Eq. p4|), the model is reduced to the 


network of identical oscillators which we discussed above. 

The approach shown above indicates there are two 
steps of the reduction scheme. The first step is the re¬ 
duction from the equation Eq. (20) to Eq. (22) or the 


continuous form Eq. (24), which means the choice of the 


Poisson manifold of the system. This is related to the the 
symmetry of the system or the recursion form of order 
parameters equations. The second step is the further op¬ 


eration from Eq. (24) to Eq. (271, which depends on the 


special choice of distribution of nonidentical parameters. 
This reduction with the Lorentzain distribution was first 
found in [5) with an ansatz for low-dimensional mani¬ 
fold, namely the OA ansatz. For the case of the Dirac’s 
delta distribution function, this reduction was discussed 
comprehensively in m with group theory analysis of the 
system. In our approach, the reduction comes from the 


same basis, i.e., the symmetry of order parameter equa¬ 
tions. 

Along the approach, we can also consider more non¬ 
identical parameters, such as the location of oscillators x 
for example in m, where the coupling function depends 
on the locations x as F(a, x), ip,u>,x, f3) with the 
distributions g(co) and q(x). The order parameter a of 
the system is defined as 




a n (cu, x)g(uj)q(x)dujdx, n = 1 , 2 ,.... (28) 


When the Fourier expansion of the coupling function for 
contains only the first three terms, for each specific oj 
and cc, the relation a n (uj,x) = (oq (oj,x)) n is exactly the 
solution for dynamical equations, which reads 


a n {u,x) = in{f 1 a n+1 (uj,x) + f 1 a n _ 1 (u},x) + f 0 al{u},x)). 

(29) 

If ct(w, x) is analytic for ui and g{u>) is the Lorentzain 
distribution, the integral Eq. (28) can be obtained for oj 
as 




a n (ui, x)g(iu)q(x)dujdx 



— i, x)q(x)dx , 


(30) 


where n = 1 , 2 ,.... 
(aq(—i, x)) n and w = 
(29), we have 


Setting 7 n (x) = a n (-i,x ) = 
—i in the dynamical equation Eq. 


li{x) =i(/i(w, 7i(x), q;)7i [x ) 2 

+ /i(w,7i (£),£*) ( 31 ) 

+ /o(w, 7 i(x),a) 7 i(a:)), 


where a = (...a_i,ao,oq...) is the order parameter of 
the system and the functional of 71 (a;) as 


ot n = J{'y\{x)) n q{x)dx, n = 1 , 2 ,.... 


(32) 


The solution of this integral differential equation Eq. (31) 


describes the structure of the local order parameter 71 ( 2 ?) 
along x. 

Up to now, we have discussed the system of all- 
connected phase oscillators in the limit of infinitely many 
oscillators and the case that only three Fourier coeffi¬ 
cients of the coupling function are nonzero. The low- 
dimensional invariant manifold, namely the Poisson man¬ 
ifold, is got for both cases. We will see in the next two 
sections that these two conditions are exactly the scope 
of the OA ansatz. Two cases beyond this cope will be 
discussed respectively. 


III. NETWORK WITH FINITE SIZE 

We have considered the general model Eq. (| 6 | in the 
limit of infinitely many oscillators, N —> 00 , in which 
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we could get the phase density of oscillators p(t, ip) and 
the corresponding continuity equation Eq. (JtJ) . Order 
parameters a n could be considered as Fourier coefficients 
of p(t, ip) , from which we get the equivalent expression 
of the dynamics of the system as the order parameter 
equations Eq. and build the approach above. 

In the limit N —> oo, the macroscopic variable, namely 
the order parameter a(t), has the limit a(t) —> a for 
steady states as the synchronous state and incoherent 
states, which gives us the basis to discuss the system an¬ 
alytically. In the case of a finite but large number of 
oscillators, i.e., N 1, the order parameter defined as 
a = (1/AO SjLi el>Pj W A1 fluctuate around the value a . 
This fluctuation depends on the number of oscillators as 
0(N~5) jTSinTI . When the fluctuation is small enough, 
as 1/V100000 ~ 0.003 for N = 100000, the collective 
behaviors of the system with finite number of oscilla¬ 
tors could be described by the approximation of infinitely 
many oscillators, for which analytical methods can be ap¬ 
plied. Some analytical methods, e.g., the self-consistent 
method and the OA ansatz, are applicable for the case 
of infinitely many oscillators or equivalently the phase 
density description. 

On the other hand, in the case of only a few number of 
oscillators, e.g., N = 10, the difference | a(t) — a | would 
be so large, as 1/vTO ~ 0.316, whose magnitude com¬ 
parable with the value of a. Obviously, it is not reliable 
to treat the system with the approximation a(t) ss a 
in this case. It is necessary to propose new approaches 
in analytically dealing with the collective behaviors of 
finite-oscillator systems. 


A. Ensemble approach 

In our approach, order parameters a n could be con¬ 
sidered as not only the Fourier coefficients of p(t, ip) but 
also the collective variables as 

1 N 

a n=N^ £intPj ' (33) 

3 =1 

which does not depend on the approximation of infinitely 
many oscillators. For the system as 

= F(a, ipj, [3),j = 1 (34) 

with the definition of order parameters Eq. ( |33| ), the 
dynamical equations for order parameters are 

K 

d « = ^E ein ^^> n = 1 ’ 2 >---- ( 35 ) 

i=i 

For the coupling function F(a, ip, (3) is 27r-periodic for <p, 
we have the Fourier expansion as 

OO 

F(a,<p,(3) = Y /,'(«,/3)e ijV - (36) 

j = ~ oo 


Substituting the expansion into Eq. (35), together with 
Eq. (33), we have the closed form of equations for order 
parameters as 


a n = m 


j = — oo 


fj n — f 5 2, . . . . 


(37) 


which is the same as we get in the limit of infinitely many 
oscillators. 

In the case of a few number of oscillators, the definition 
of order parameters can also be considered as the coor¬ 
dinate transformation that transforms the microscopic 
variables {<Pj} to macroscopic variables {a n }. Therefore 
it is inspiring that the dynamics of phase oscillators with 
corresponding initial conditions {</?., (0)} is equivalent to 
the dynamics of order parameters with corresponding ini¬ 
tial conditions {a n (0)}. The number of oscillators, no 
matter finite or infinite, has no influence on the dynam¬ 
ical equations of order parameters. This forms the im¬ 
portant basis of our approach. 

Following this approach, when only the first three 
Fourier coefficients of the coupling function are nonzero, 
by setting a n = a™ and substituting the relation to the 
dynamical equations, the equations Eq. (371 for all the 
a n will be reduced to a single one as 


di — i(/i(aq)aq + /i(aq) + /o(aq)aq). (38) 


It appears that we could get the Poisson manifold again 
even the size of the system is finite. However, the finite- 
size effect should be taken into account with care. 

As a matter of fact, in the case of a finite number of 
oscillators, the Poisson manifold with a n = a" is not 
attainable for the system. To see this, take the system of 
N = 2 as an example. For the first two order parameters 
aq and a. 2 , by definition 


Oil 



j =i 


a.2 


1 

2 


2 


Y 0 2iVj 


(39) 


if the system described by oq, q ?2 is in the Poisson mani¬ 
fold, the relation a 2 = af could be rewritten in terms of 
the definition Eq. (f33l) as 


\Y e2iVi = (bE e ^) 2 - 


(40) 


1=1 


3 =1 


However, the equality Eq. (40) is not naturally valid. 


By using a simple calculation of the difference between 
the left-hand and the right-hand terms in Eq. (40) one 
obtains 

2 2 

bE e2 ^ - (^E ^) 2 = j(e^-e^) 2 . (41) 


3 = 1 


3 =1 


This indicates that the system can evolve on the Poisson 
manifold only when ip 2 = ipi■ In general, the relation 









6 


a n = Oil gives exactly infinitely many independent con¬ 
straints like Eq. (40) with n = 2,3,..., and any solutions 


for finite oscillators will be determined as the trivial one 
as (fj = tpi for all the index j, which means that ex¬ 
cept the synchronous state, all the states of the system 
of finite oscillators lie out of this Poisson manifold. The 
system can only evolve on the Poisson manifold in the 
limit of infinitely many oscillators. 


Whereas, as a matter of fact, the equations Eqs. (37) 


hold for all N, whether the size of system is finite or 
infinite, and the synchronous state always satisfies the 
relation of the manifold. The Poisson manifold can be 
used at least approximately in the vicinity of synchronous 
state. However, we need some new methods and approx¬ 
imations. 

Let us consider an ensemble with M identical systems 
of N oscillators which have the same dynamical equations 
and same parameters. The initial phases of oscillators are 
chosen from the same distribution, which makes the mean 
values of order parameters a n over the ensemble have the 
limit when M —> oo. This leads to the definition of the 
ensemble order parameter (a n ) as 


1 


M 


(d n ) = Jim Yr^oi^\n = 1,2,..., (42) 

M—>oo M z ' 

3 =1 

where M is the number of sampling systems in the en¬ 
semble and is the nth order parameter for the jth 
system in the ensemble. Taking the ensemble average for 
both sides of the dynamical equations Eq. (37), we get 


(a n ) = m((/ia n+ i) + (/ia n _i) + (foa n )), (43) 

where n = 1 , 2 ,..., and (•) means ensemble average. 

In general Eq. (43) is not solvable because the terms 
in the right side as (/ia n +i) cannot be simply described 
by (a n ) in general. However, if fj is independent of all 
the systems in the ensemble as (fjOc n ) = fj{a n ), then we 
can get exactly the dynamical equations of the ensemble 
order parameters, as 


d{oi n ) 

dt 


= m(/i(a„+i) + fi(a n -i) + f 0 (a n )), (44) 


where n = 1,2 ,... and the terms /±i,o may depend on 
(a n ). Similar to Eq. (37), there is a solution for the 
ensemble order parameter equations Eqs. (44) as (a n ) = 
(ai) ra , which is exactly the Poisson manifold. 

Moreover, for more general cases ( fja n ) ^ fj{a n ), but 
when the terms (fjOt n ) can be approximated by ( fj){ot n ), 
namely the statistical independence, we can also get the 
dynamical equations of the ensemble order parameters 
(a n ) with similar approach. The validity of this approach 
can be measured by the error terms 


e jn — ( fj a n ) ( fj)( a n ), 


(45) 


B. The star Sakaguchi-Kuramoto model 

In the following, let us take the star Sakaguchi- 
Kuramoto model as an example, which is a typical topol¬ 
ogy and model in grasping the essential properties of het¬ 
erogeneous networks and synchronization process, as 

N 

0 h = U h + A ^2 sin(0j - 0 h - /?), 

j =i ( 46 ) 

0j = uj + X sin (Oh — 0j — p), 

where 1 < j < N, uih,0h and u>j,0j are the natural fre¬ 
quency and the phase of hub and leaf nodes respectively, 
A is the coupling strength and p is the phase shift. By 
introducing the phase difference <pj = Oh — 0j, the dy¬ 
namical equation can be transformed into 

N 

<Pj = Ac o — X ^ sin(y>fc + P) — X sin(^- — P), (47) 

k =1 

where A to = tOh — ojj is the natural-frequency difference 
between hub and leaf nodes. 

By introducing the order parameter re 1 * = a = 
jj eliPj j w e could rewrite the dynamics as 

l p j = f e m +g + f e -m : j = (48) 

where / = i^e -1/3 , g = Aw — AiVr sin(4) + p). Hence the 
system is in the framework which we discussed above, 
with the first three nonzero Fourier coefficients of the 
coupling function. With our approach, the dynamical 
equations for order parameters can be obtained as 

din = in{fa n+ i + fa n -1 + ga n ),n = 1,2,- (49) 

In this specific model, we cannot simply make the ap¬ 
proximation as N —> oo because g = Aw — AlVr sin($+/3) 
diverges with N —» oo. There are two ways in further 
simplifying the system. First, we could set a hypotheti¬ 
cal system as 


M 

flj = Aw — NX ^ 
j=i 


1 

M 


sin(i?j + P) 


X sin(i?j 


P), (50) 


where 1 < j < M, dj is the hypothetical phase oscil¬ 
lator and N is considered as a parameter in this sys¬ 
tem. Defining the order parameter for this system as 
re 1 * = a = YlfjLi ell9j > Eq. (50) becomes 


<i>3 = f e ' Vj + 9 + fe lLpj , j = !,■■■ ,M, 


(51) 


where / = i^e _1/3 , g = Aw — AA r rsin(4> + /?). Following 
our approach, the corresponding order parameter equa¬ 
tions read 


din — k^(/' Oinj-1 T f Oi n — 1 -\~ gOin'), Tl — 1,2,..., 


with j = ±1,0 and n = 1,2 ,.... 


(52) 
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which are exactly the same as Eq. (49). But for this 


hypothetical system, we could take the limit of infinitely 
many oscillators as M —> oo, which could be discussed 
analytically with the traditional OA anstaz. 

Obviously, these two systems are quite different for 
phase variables, one with finite oscillators the other with 
infinitely many oscillators, but they have the same order 
parameter equations. In the following we will see that the 
hypothetical system Eq. (50) is exactly a representation 


of the ensemble of the original model Eq. (47). 

On the other hand, following our ensemble approach, 
for the system Eq. (47), let us choose an ensemble con¬ 


sisting of systems with the same parameters A, A u and 
N , and different initial conditions chosen from the same 
distribution, as the Poisson kernel. We have the ensemble 
average of the dynamical equation as 

(d„) = in(i^e _1/3 (a n+ i) - i^e I/3 (a„_i) + Aw(ct n ) 


+ -AA((a 1 a n )e I/3 “ («i a n)e 1/3 )), 

where n = 1,2,.... Setting 

(oqcr n ) — (oq)(cr n ) -f- Cni, 
{a\a n ) = ( a\){a n ) + e n 2, 


(53) 


(54) 




FIG. 1: (color online) (a) and (b) The first three error terms 
E n and E n versus time for the system with N = 10. (c) and 
(d) Order parameters getting from numerical simulation and 
approximated OA ansatz. Results from approximated OA 
ansatz is the red line, and every single simulation results is 
the blue lines in (c) with ensemble average the light blue line 
in (d). 


for n = 1,2,.... If \e n i\ < |(ai)(a n )| and |e n2 | < 
|(di)(a„)|, or typically \e n i\ < 1 and |e„ 2 | < 1, they 
can be seen as the perturbation terms for the dynamics 
of ensemble order parameters in Eq. (53). Ignoring the 
perturbations, we could get 


i) = in(i^e lp (a n+ i) - i^e I/3 (a„_i) 
A u{a n ) + ^A N{{a 1 ){a n )e' ft 


(55) 


-(arXa^e-*)), 


and this is the same as Eq. ( [52] ) for hypothetical system. 
In this case, the ensemble order parameters are the same 
as the order parameters for the hypothetical system. 

For the case of a finite but large number of oscillators, 
the hypothetical system is exactly the continuous model 
which shares the same order parameter equations but 
with infinitely many oscillators. As in the case of a few 
oscillators, the hypothetical system is a presentation of 
the ensemble for the systems, where the difference terms 
e„i and e n2 by this method will not only introduce some 
fluctuations but also some systematic errors, which can 
be described by 


E n — |(oqo; n ) (ou)(o; n )|, 

, _ (56) 

E n — |(oqQ! n ) (ou)(o; n )|, 

with n = 1,2 ,... for this specific model. 

In terms of numerical computation, we can check the 
above approximation by examining E n and E n for the 
ensemble of M systems. The approximation E n and 


E n can be divided into two parts as E n = /„ + s n and 
E n = f n + s n where f n and f n are the fluctuation parts 
proportional to 1 /\fM which depend on the size of the 
ensemble and s n and s n are the systematic parts which 
are introduced by the statistical independence assump¬ 
tion. 

Take N = 10 as an example. When M = 1000 the 
fluctuation parts are ignorable. We find that except for 
E 1 is around 0.03, all of E n and E n for n < 3 are much 
smaller than 0.01 in this case as plotted in Fig. 0 a) and 
(b). Comparing the typical magnitude of values of order 
parameters as 10 1 , the approximation of the ensemble 
approach is obviously reasonable. 

For the initial conditions chosen from the Poisson ker¬ 


nel, with the dynamical equations Eq. (55), the ensemble 


order parameters will evolve on an invariant manifold, 
i.e., the Poisson manifold with ( a n ) = ( a\) n . Denote 
(ai) by z, we have 




+ i(Aw + N(ze i/3 - ze~ i0 ))z + 




( 5? ) 

The difference between z and the ensemble average of nu¬ 
merical simulations is checked in Fig. m , which shows 
that it is reasonable to describe the behavior of the en¬ 
semble by Moreover, for every single system in the en¬ 
semble, with some fluctuation, it can also be described by 
the ensemble order parameter and hence by z, as shown 
in Fig. Qc). 

With this ensemble assumption works, the two- 
dimensional equation Eq. (57) in the bounded space 
as \z\ < 1 describes the dynamics of the coupled os- 





































cillator system. Every stationary state for phase os¬ 
cillators system has its counterpart in this space for 
z = re 1$ . For example, the in-phase state(IPS) de¬ 
fined as <Pj(t) = <p q (t) = (ft corresponds to a limit cy¬ 
cle as r(t) = 1. The synchronous state (SS) defined as 
< fj(t ) = (f q (t) = (f corresponds to a fixed point with 
r(t) = 1 and $(i) = c. The splay state(SPS) defined 
as (fj{t ) = ip(t + with T the period of <p(t) corre¬ 
sponds to a fixed point with r(t) = ro < 1 and $(f) = c. 
In the two-dimensional order-parameter space, it is easy 
to analyze the existence and stability conditions of these 
states, and the basion of attraction for each state in the 
coexistence region can also be conveniently determined 
|19| . This is the reason why we introduce the invariant 
manifold governed by low-dimensional dynamics. 

As we see, in the case with only finite oscillators, the 
ensemble order parameter plays the dominant role in re¬ 
vealing the key collective structure of the system. The 
price of the ensemble order parameter description is the 
error induced by the statistical independence assump¬ 
tion, whose validity can be checked by numerical results. 
In the next section we will see that when more than three 
Fourier coefficients of the coupling function are nonzero, 
where the traditional order parameter scheme fails, one 
can still discuss the collective dynamics in terms of the 
ensemble order parameter approach. 


IV. COMPLEX COUPLING FUNCTION 


as the transformation of Eq. (581 with ip = m4> and 
fj = gj.m■ The order parameters a for tp is related to 
the order parameters a' for <j> as 


K K 

— _ \ '' piruflj _ _ \ " _ r 

zy u-n m- 


(62) 


i=i j=1 

And the phase density of <fr with the transformation reads 

= p(v) 


= - E 


= — 'V a 1 e y 

2^ a 3-m e 


2tt 


i j-m(f> 


(63) 


27T 


This indicates that all the terms a' n with \n\ ^ j ■ m,j = 
0,1,2... are zero. Moreover, we can transform the Pois¬ 
son manifold a n = (oq) n for <p to the invariant manifold 
a n-m = ( a m) n f° r </*) which is governed by 


“m — iTO (/m( a m) 2 + fm + /S a m)i 


(64) 


corresponding to Eq. (601 for aq. This is exactly the low¬ 


dimensional manifold for the coupling function Eq. (61) 


with higher order Fourier coefficients. This manifold is 
based on the mth order parameter a' m and the statement 
that all the terms a' n with \n\ ^ j • m,j = 0, 1 , 2 ... are 
zero. This is a bit queer in this respect. 

Let us take the case m = 2 as an example [T5] . In this 
case only f ±2 and fo among the Fourier components of 
the coupling function are nonzero, i.e., 


In the above sections, for the coupling functions with 
only first three terms as 

F(a,(p) = fi(a)e llp + /_i(a)e _1V3 + f 0 (a), (58) 

we have the dynamical equations for order parameters 
a n as 

d n = m(/ia„+i + fia n -i + f 0 a n ),n = 1,2,- (59) 

On the invariant manifold a n = a™, all these equations 
are reduced to a single one 

oq = i(/i(oq)aq + /i(aq) + /o(aq)aq). (60) 

For general situations, the coupling function can not be 
simply truncated to only the first terms. It is thus impor¬ 
tant to consider the order parameter approach for more 
complicated cases. 


A. Higher order coupling 


Firstly, let us consider a coupling function with higher 
order Fourier coefficients like 

G(c/, </>) = g m (ot')e irn ^ + g„ m ( a ')e- im * + g 0 (a'), (61) 

where to > 1 is an integer and a' are the order parame¬ 
ters for the system. Note that Eq. (61) can be regarded 


F(a,p) = f 2 (a)e 2ltp + f_ 2 (a)e 2lip + f 0 (a). (65) 

According to the above analysis, the system has the in¬ 
variant manifold as cc2n+i = 0,O!2n+2 = a 2 +1 > n — 0, or 
equivalently a n = l/2aE[(—l) n + 1], with dynamics 

«2 = 2i(/ 2 «2 + h + /oM- (66) 

The phase density of oscillators in this case is 

, OO 

P(<P) = ^ E (67) 

j= — c» 

By using our order-parameter-analysis approach, for 
the case of to = 2, it is not hard to get the the equations 
of order parameters a n as 

d2ri+l = i(2n + l)(/2d;2n+3 + /2C*2n-l + /oO ? 2n+l) ) 

Ct-ln = i(2n)(/2Qf2n+2 + /2Ct2ra-2 + foOt2n), 

(68) 

with n = 1,2,.... The manifold with oi2rt+i = 0, a2n = 
a .2 ,n > 0 is obviously the solution of these equations. 
This invariant manifold separates the relation a n = a" 
into two parts, where the part 02 n = ct 2 conserves the re¬ 
lation, and the odd part Q!2n+i = 0 can be regarded as a 
special choice of the relation. Even if in the state defined 
as \ot 2 n | = 1021 n = 1, we still have aq = 0, which repre¬ 
sents the state of cluster synchrony as discussed in T81 . 
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And the manifold with a 2 n+i = 0 ,tt 2 n = ot 2 ,n > 0 can 
be regarded as a sub-manifold of the Poisson manifold. 

An interesting issue here is the reason why we set all 
the a 2 n+i = 0 other than a 2 n+i = ay" -1 " 1 . To see this, 


we just substitute the relation « 2 n+i = a 2n+1 ,n > 0 


into the odd part of Eq. ( 68 ), then the equations will be 


reduced to the following equations, 

£*1 = i(Mi + Mi + fo&i), 
on = i(h a i + h a i 1 + fo&i), 


(69) 


where the first equation is got from n = 0 in Eq. ( 68 ), 


and the second equation is got from all the others n > 0 


in Eq. ( 68 ). Comparing the two equations in Eq. (69), 


the consistent condition requires dq = , and this im¬ 

plies that |oq| = 1. Hence, there isn’t a general relation 
of ct 2 n+i, apart from the vicinity of synchronous state. 
Thus we should consider another special solution for the 
odd part of Eq. ( 68 ) as a^n+i = 0, n > 0. 


On the other hand, for the even part of Eq. ( 68 ), by 


setting a 2 n = Q%, the equations will be reduced to a 
single one as 


d 2 = 2i(f 2 al + f 2 + /oa 2 ). 


(70) 


The conjugate part a _2 is separated from 02 naturally in 
this equation, and we can get the low-dimensional mani¬ 
fold governed by the order parameter equation Eq. (70), 
as a sub-manifold of the Poisson manifold. 


B. The Dominating-term assumption 

As mentioned above, the Poisson manifold no longer 
exists when the order parameter relation a n = a” can 
not be used to reduce the order parameter equations to 
the single one. For the case of only higher terms this 
problem can be fixed by separating the order parameters 
into groups and setting some of them zero with which a 
special sub-manifold of the Poisson manifold can be used 
to get low-dimensional behaviors. However, when more 
than three terms in Fourier expansion are considered, 
even the sub-manifold on longer exists. 

Consider a coupling function with the first five nonzero 
Fourier coefficients, as 


F(a,'p)= ^ f n (a)e mv , 


(71) 


the corresponding equations of order parameters read 
now 

a n = m(/ 2 <a :„ + 2 + / 2 a n -2 + fi&n+i + fiOi n -i + /o Q n)i 

(72) 

with a- n = a n , n = 1,2,.... In this case, if we suppose 
a n = a", then two different dynamical equations will be 
obtained 


Oi 1 — i(Mi + f 2 oi\ + /ioy + fi + /oaq), 
on = i(/2Q^i + Mi 1 + /i«i + fi + fo a i)- 


(73) 


Hence we will have a_i = aj” as the requirement of 
coincidence of these two equations. The relation of a n = 
a™,n > 0 is broken by the conjugate terms and there is 
no way in finding out a separated solution as we did for 
Eq. ( 681 . The only solutions with the relation a n = a" 
for this system seems to be either a„ = a" = 0 or \a n \ = 

|ai|" = 1. None of them are expected for our analysis, 
unless in the vicinity of synchronous state with |oq| = 1 
or in the vicinity of the incoherent state with oq = 0. 

On the other hand, for n > 2 all the dynamical equa¬ 


tions for the nth order parameter in Eq. (72) are reduced 
to the same equation by a n = a", and the symmetry of 
Eq. (72) can be represented by a n = a™ at least for all 
the order parameter equations for n > 2. Moreover, due 
to the bound |aq| < 1 , the higher the order components 
a n become smaller with increasing n. If we choose the 
initial conditions as a n = a™, then along the evolution 
of Eq. (72), the differences of higher terms a n with a™ 


will keep small enough. By ignoring the differences and 


making the approximation that a n 
the dynamical equations for oq as 


a7, we can obtain 


Oil — i(/ 2 «i + / 2 od + /ioq + fi + /ocki), 


(74) 


which can be regarded as the main term of the dynamics 
Eq. (72). From now on we call this approach dominating- 


term assumption. We will further show that this approx¬ 
imation is pretty efficient in dealing with systems with 
complicated coupling functions. 


Let us take the system described by Eq. (47) as an 


example. By considering higher terms, we have the dy¬ 
namical equations 

<Pj = K/ie^ + g + fie-w + f 2 e 2 ^ + f 2 e~ 2i (75) 

where j = 1, • • • , N and jq = i4^, f 2 = i4f, g = 
Aw — ATVrsin(<I>) with re 1$ = aq. By introducing the 
dominating-term assumption, the dynamical equation for 
oq reads 

01 = i{hoi\ + Mi + /M + /1 + gai). (76) 


If the dominating-term assumption works, Eq. (76) can 


be used to study the collective behaviors of the coupled 
oscillators, which can be regarded as a sort of approxi¬ 
mate low-dimensional behaviors. 

The validity of the dominating-term assumption can be 
checked via numerical simulations. For the initial state 
chosen from the Poisson kernel distribution, by using nu¬ 
merical simulation, we have the order parameter from the 
continuity equation of the phase density and the approx¬ 
imate manifold. In Fig. [2ja) the order parameter oq is 
plotted. It can be found that the result in terms of the 
dominating-term approximation coincides very well with 
numerical results. In Fig. |2](b), the density of oscillators 
at a given time is also plotted, and the result shows that 
though the distribution is indeed no longer the Poisson 
kernel distribution, it is instead close to it. This gives us 
the mechanism of how this assumption works. 
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FIG. 2: (color online)(a) Order parameters getting from nu¬ 
merical simulation (red line) and approximate OA ansatz 
Eq.( |76| l(blue line), (b) Distribution of phase variables getting 
form numerical simulation (red) and reconstituting Poisson 
kernel with order parameter ai (blue). The system consid¬ 
ered as A = Ai = 0.4, A 2 = 0.3, N = 1, Acu = 5 


According to numerical results, we can see that the 
dominating-term assumption works pretty well, and the 
system shows approximated low-dimensional behaviors. 
With the dynamical equations for the first order param¬ 
eter Eq. (76), we can do the study further to find the 
property of collective behaviors of the system, which will 
be discussed in detail in another paper. 


V. DISCUSSION 

In this paper, we focused on the theoretical description 
of low-dimensional order-parameter dynamics of the col¬ 
lective behaviors in coupled phase oscillators. We have 
derived the closed form of dynamical equations for or¬ 
der parameters, from which we find the Poisson kernel 
as an invariant manifold and the well-known OA ansatz. 
Different from the traditional Ott-Antonsen ansatz and 
Watanabe-Strogatz’s approach, our approach is suitable 
for systems of both finite and infinite, identical and non¬ 
identical oscillators. In our approach, the scope of the 
OA ansatz is determined as two parts, i.e., the limit of 
infinitely many oscillators and the condition that only the 
first three Fourier coefficients of the coupling strength are 
nonzero. 

By using our order parameter analysis, we also dis¬ 
cussed two cases that go beyond the scope of the OA 
ansatz, i.e., the case of a finite-size system and the case 
of coupled oscillator systems with more complicated cou¬ 
pling functions. We have discussed the reasons why the 
OA ansatz cannot be used directly to these two cases 
and we further developed the approximation methods 
to deal with these difficulties with the OA ansatz. We 


developed two methods, namely the ensemble method 
and dominating-term assumption. It is shown that these 
schemes work pretty well, and their validity has been 
checked by the numerical simulations. 

For the case of a finite number of coupled oscillators, 
it is shown that the system is out of the domain of 
low-dimensional manifold of the coupled order parameter 
equations. Hence, from the view of the geometry theory 
of manifolds, the finite size of system brings some fluctu¬ 
ations along the manifold which is introduced from initial 
states. If the dynamics of order parameter equations is 
not too complicated and the manifold satisfying a n = a" 
can be used to describe the dynamics approximately, all 
the trajectories from different initial conditions around 
the Poisson manifold will be described by the mechanism 
of this manifold approximately. This is the meaning of 
the ensemble method and the reason for why it works in 
some models. 

In other case, with a more complicated coupling func¬ 
tion, the Poisson manifold is not a solution of the dy¬ 
namical system. However, for the infinitely many cou¬ 
pled equations for a„, except for few of them, the re¬ 
lation a n = a” indeed reflects the symmetry of these 
order parameter equations. When we choose the initial 
states in this manifold, the evolution of the dynamics will 
certainly be influenced by the symmetry of it, which is 
described by a n = a" approximately. The method with 
dominating-term assumption is indeed designed to show 
this effect, and consequently the system is simplified by 
it for some models. 

In this paper, we have developed the order param¬ 
eter analysis, with which we get the low-dimensional 
behaviors of coupled phase oscillators, such as the OA 
ansatz, the ensemble order parameter approach and the 
dominating-term approximation. However, to fully un¬ 
derstand the low-dimensional collective behaviors we 
need a more comprehensive understanding of the Pois¬ 
son manifold and its relation with the dynamics of the 
system, which is still an open question. We believe that 
the order parameter analysis should be a powerful tool 
in helping to reveal the mechanism of low-dimensional 
collective behaviors. 
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